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Recent  developments  In  computational  fluid  mechanics  and  grid  generation 
have  made  It  possible  to  consider  complex  three-dimensional  laminar  and  turbu¬ 
lent  flows  using  nonorthogonal  coordinates.  However,  the  necessary  equations 
In  generalized  nonorthogonal  coordinates  are  not  readily  available  In  a  single 
reference.  The  purpose  of  this  report  Is  to  present  these  equations  for  an 
incompressible  fluid.  Methods  of  curvilinear  tensor  analysis  are  used  to 
develop  the  most  commonly  used  equations.  Their  relationship  with  the  more 
familiar  equations  in  orthogonal  coordinates  is  shown.  Finally,  the  practical 
problems  of  evaluating  the  many  geometric  coefficients  which  appear  in  these 
equations  are  examined  by  considering  numerically-generated  coordinates  for 
flow  around  three-dimensional  bodies. 
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turbulence  model  constants 
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I.  INTRODUCTION 


A  problem  that  frequently  arises  In  the  development  and  application  of 
general  numerical  methods  for  the  solution  of  fluid-flow  equations  Is  the 
choice  of  a  coordinate  system.  This  problem  becomes  particularly  crucial  In 
the  analysis  of  three-dimensional  flows,  whether  Internal  or  external,  since 
the  coordinates  chosen  must  satisfy  many  conflicting  requirements  and  influ¬ 
ence  the  accuracy  of  the  solutions.  For  example,  the  application  of  boundary 
conditions  on  a  solid  wall  Is  greatly  facilitated  by  locally-orthogonal  body- 
fitted  coordinates,  while  proper  resolution  of  the  flow  In  certain  regions 
necessitates  the  use  of  nonorthogonal  curvilinear  coordinates. 

Many  methods  are  now  available  for  the  construction  of  general  nonortho¬ 
gonal  coordinates,  but  their  use  In  the  actual  solutions  of  fluid-flow  prob¬ 
lems  Involving  three-dimensional  geometries  is  still  in  Its  Infancy.  Once  a 
coordinate  system  Is  selected  for  a  given  geometry,  or  class  of  problems, 
there  still  remains  the  task  of  transforming  the  equations  of  motion  from  the 
familiar  orthogonal  coordinates  to  those  In  the  chosen  system.  Two  different 
approaches  have  been  adopted  for  this  purpose.  One  of  these  uses  partial 
transformations,  In  which  only  the  Independent  coordinate  variables  are  trans¬ 
formed,  leaving  the  dependent  variables  (l.e.  the  velocity  components)  In  the 
original  orthogonal  coordinates.  This  approach,  which  has  been  used  by  Chen 
and  Patel  (1985a,  1985b),  among  others,  has  the  advantages  that  the  resulting 
equations  are  relatively  simple  and  the  results  can  be  readily  Interpreted. 
If  the  physical  geometry  Is  highly  curved,  such  an  approach  may  lead  to 
increased  numerical  diffusion  due  to  large  skew  angles  between  the  velocity 
components  and  the  faces  of  the  computational  cell.  The  alternative  Is  to 
transform  the  equations  completely.  Including  the  dependent  as  well  as  the 
independent  variables.  This  approach  has  been  used  principally  In  two- 
dimensional  flow  problems  since,  In  three  dimensions,  the  equations  involve 
many  more  geometric  coefficients.  This  leads  to  Increased  storage  require¬ 
ments  if  the  coefficients  are  calculated  once  and  for  all,  or  longer  execution 
times  If  they  are  calculated  each  time  they  are  needed.  It  is  also  possible 
that  numerical  errors  associated  with  the  evaluation  of  the  geometric  coeffi¬ 
cients  could  adversely  affect  the  solution  of  the  flow  equations.  However, 


use  of  a  complete  transformation  provides  a  method  to  obtain  solutions  for  a 
much  wider  variety  of  problems  than  Is  possible  with  the  partially-transformed 
equations. 

In  spite  of  the  rapid  developments  of  methods  for  the  construction  of 
generalized,  nonorthogonal ,  curvilinear  coordinates  for  fluid  flows  with 
three-dimensional  geometries  (see,  for  example,  Thompson,  1982),  there  does 
not  appear  to  be  a  ready  reference  in  which  one  can  find  the  corresponding 
equations  for  laminar  and  turbulent  flows.  In  most  Instances,  the  equations 
are  derived  as  needed  for  a  particular  coordinate  system  and  it  is  usually 
quite  difficult  to  make  generalizations  or  extensions. 

The  partial  transformation  of  the  Reynolds-averaged  Navler-Stokes  equa¬ 
tions,  and  the  transport  equations  of  the  well  known  k-e  model  of  turbulence, 
from  a  general  curvilinear  orthogonal  coordinate  system  to  a  nonorthogonal 
system  were  presented  by  Chen  and  Patel  (1985a).  The  main  purpose  of  the 
present  report  Is  to  document,  for  general  reference,  the  fully-transformed 
equations.  Because  methods  of  tensor  analysis  are  used  throughout  the  report, 
a  summary  of  the  basic  concepts  of  tensor  analysis  is  presented  first.  The 
Reynolds-averaged  Navler-Stokes  equations  and  the  transport  equations  for  the 
Reynolds  stresses  are  then  developed  for  an  Incompressible  fluid  without 
reference  to  any  particular  turbulence  model.  The  equations  of  turbulent 
kinetic  energy  and  Isotropic  dissipation,  which  form  the  basis  of  the  k-e 
model,  are  then  presented.  The  Interrelationship  between  coordinates  is  then 
demonstrated  by  reducing  the  equations  In  generalized  coordinates  to  those 
given  by  Nash  and  Patel  (1972)  for  orthogonal  coordinates.  Finally,  to 
provide  practical  examples,  typical  geometric  coefficients  that  appear  In  the 
equations  are  evaluated  from  numerically-generated  nonorthogonal  coordinates 
for  two  shlp-llke  three-dimensional  bodies. 


II.  ELEMENTS  OF  TENSOK  ANALYSIS 


The  tools  of  curvilinear  tensor  analysis  will  be  used  extensively 
throughout  this  report  because  they  provide  the  most  direct  method  by  which 
equations  valid  for  nonorthogonal  coordinate  systems  may  be  derived.  The 
texts  by  Lass  (1950)  and  Arls  (1962)  are  the  main  references  for  the  material 
concerning  curvilinear  tensors  presented  herein.  This  section  provides  a  sum¬ 
mary  of  useful  concepts  and  expressions  frequently  used  In  the  later  sections 
and  which  the  reader  may  find  useful  In  the  application  or  extension  of  the 
results  presented  In  this  report. 

II. 1.  Coordinate  Transforaatlons 

Curvilinear  tensor  analysis  differs  fundamentally  from  Cartesian  tensor 

analysis  In  that  the  base  vectors  (unit  vectors)  can  be  functions  of  position. 

1  12  3 

General  curvilinear  coordinates  are  denoted  by  x  =  (x  ,  x  ,  x  )  and 

1  12  3 

5  =  (C  ,  C  ,  ).  The  coordinates  can  be  defined  either  analytically  or 
numerically  and  may  be,  In  general,  orthogonal  or  nonorthogonal.  The  follow¬ 
ing  summation  convention  Is  used  throughout  this  report;  If  an  Index  appears 
In  a  product  of  terms  once  as  a  subscript  and  once  as  a  superscript  a 
summation  upon  the  index  from  1  to  n,  where  n  is  the  dimension  of  the  space. 
Is  Implied.  For  example.  In  a  three-dimensional  space 

112  3 

A  Aj  =  A  Aj  +  A  A2  +  A 

where  1  Is  a  dummy  index.  A  free  Index  Is  an  Index  which  is  not  repeated  In 
the  above  manner.  In  order  to  avoid  confusion,  the  same  letter  should  not  be 
used  to  denote  more  than  one  free  Index  In  an  equation.  The  consistency  of  an 
equation  may  be  quickly  checked  by  noting  whether  the  free  Indices  appear  in 
the  same  position  (superscript  or  subscript)  throughout  the  equation.  For 
instance. 


A^B^  = 


E 

k 


Is  consistent,  hut 


-3- 


Is  not  consistent 


A  general  transformation  relating  two  coordinate  systems  o^,  Is 


defined  by 


,  ,  1  2  nv 

f  (a  ,  o  ), 


1  *  ly2|a«»fn 


(2-1) 


where  f  could  either  be  an  analytic  or  numerical  transformation.  In  order  to 
ensure  that  the  transformation  be  a  proper  one-to-one  Invertible  mapping  we 
should  be  able  to  solve  for  the  such  that 

=  f“^  (B^  B^,...,b”),  1  -  1,2, ...,n  (2-2) 


The  conditions  under  which  f~^  exists  are  determined  by  the  Implicit  function 
theorem  (Bartle,  1975;  Hildebrand,  1976).  A  transformation  Is  said  to  be 
proper  if  the  Jacobian  determinant  does  not  vanish.  For  a  three-dimensional 
space  J  is  of  the  form 


36^ 

3B^ 

3B^ 

1 

2 

3 

3a 

3a 

3a 

2 

2 

2 

3B 

3B 

3B 

^  1 

,  2 

—J 

3a 

3a 

3a 

3 

3 

3B 
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3B 

1 

2 

3 

3a 

3a 

3a 

Recall  from  calculus  that  the  Jacobian  represents  the  ratio  of  elemental 
volumes  for  a  change  of  variables  In  a  volume  Integral,  l.e. 


dV 


3(8^  B^ . a")  .  , 

1  2  n, 

9  (oi  y  0  I  •  •  •  jOt  ) 


(2-4) 


where  dV^  is  the  Initial  elementary  volume.  Should  the  Jacobian  of  the  trans¬ 
formation  vanish  an  elemental  volume  would  be  mapped  Into  a  point  or  vice 
versa.  Clearly  this  violates  the  requirement  that  the  transformation  be  a 
one-to-one  mapping;  thus  we  should  consider  only  proper  transformations. 


Il»2.  Covariant  and  Contravarlant  Vectors 


Consider  a  proper  coordinate  transformation  of  the  form 

1  ■  1,2,. ..,n  (2-5) 


Now  the  differential  along  a  coordinate  direction  Is  given  by 


dr- 


^dxJ 


(2-6) 


ax-' 


where  a  summation  on  J  Is  Implied.  A  vector  whose  components  transform  In 
this  manner  Is  said  to  be  a  contravarlant  vector.  The  general  transformation 
rule  for  contravarlant  vectors  Is  given  by 


A 


1 


ax^ 


(2-7) 


Differentiating  the  transformation  of  the  coordinate  differential  (2-6)  with 
respect  to  time  yields 


dii 

dt 


a?*  dx^ 


(2-8) 


i 

where  jp-  »  V-'  are  the  components  of  the  velocity  vector.  Thus,  the  velocity 
Is  a  contravarlant  vector,  as  are  the  acceleration  and  force  vectors. 
Consider  a  scalar  function 


.X  1  2  n 

<Ji(x  ,  x  ,...  ,x  ) 


whose  partial  derivative  with  respect  to  a  new  coordinate  system  is 


^  3<t>  ax^ 
a?^  ax^  ac^ 


(2-9) 


A  vector  whose  components  transform  In  this  way  Is  said  to  be  a  covariant 

vector.  Observe  that  are  the  components  of  the  gradlant  vector.  The 

ax^ 

general  transformation  rule  for  covariant  vectors  Is 


U-IOJ 


The  concept  of  contravarlant  and  covariant  vectors  can  also  be  Introduced 
by  considering  the  differential  of  a  position  vector.  Any  line  element  vector 
dr  Is  related  to  the  differential  change  of  coordinate  d?^  hy 


(2-11) 


where  r  Is  the  position  vector,  and  the  expression  Is  summed  over  the  values 

3r 

of  1.  It  Is  convenient  to  define  the  derivative  — ^  ss  the  ith  covariant  base 

35 

vector  a^  ,  l.e. 


(2-12) 


so  that  the  line  element  vector  can  be  written  as 


dr  ••  d5  a, 
~  ~1 


(2-13) 


More  generally,  any  vector  A  can  be  decomposed  In  the  same  manner  as 


A  =  A  a. 


(2-14) 


where  A^  denotes  the  contravarlant  components  of  the  vector  A.  The  reciprocal 
(contravarlant)  base  vectors  are  defined  by 


5  =  j  (Sj 


(2-15) 


with  (l,j,k)  In  cyclic  order  (e.g,.  If  1  «  1  then  j  •  2  and  k  =  3) ,  such  that 


i  A 
•  ■  i,  ■ 


(2-16) 


where  the  Kronecker  delta  6,  la  a  mixed  second-rank  tensor  defined  by 


and  J  is  the  Jacobian  defined  by 


J  -  a  •  (a  X  a. )  (2-18) 

~1  ~j  ~k. 

For  convenience,  we  shall  define 

b  »  a,  X  a,  *  Ja 

~J  ~k 

So  that  the  1th  component  of  the  1th  contravarlant  base  vector  can  be 
written  as 


.1  m  n 

b.  =  e.  a,  a, 
I  Imn  j  k 


(2-19) 


where  (1 


,j,k)  Is  In  cyclic  order  and  the  alternating  tensor  e  Is 

Imn 


£ 


Imn 


'  I  If  Imn  Is  an  even  permutation  of  123 
<  -1  if  Imn  Is  an  odd  permutation  of  123 
^  0  If  any  two  indices  are  equal 


(2-20) 


A  vector  A  could  also  be  decomposed  In  terms  of  the  contravarlant  base  vectors 
1 

a  as 


A  -  A.a^  (2-21) 

where  Aj^  are  the  covariant  components  of  the  vector  A.  The  covariant  and 
contravarlant  base  vectors  are  not  usually  of  unit  length,  but  a  set  of  unit 
vectors  tangent  and  normal  to  the  coordinate  lines  can  be  defined  as 

(2-22a) 

(2-22b) 


respectively^  The  contravarlant  components  of  a  vector,  for  example  velocity, 
point  along  the  tangent  to  coordinate  curves.  Covariant  components,  such  as 
the  gradient  of  a  scalar,  point  in  the  direction  normal  to  coordinate  sur¬ 
faces.  The  orientation  of  the  base  vectors  are  shown  in  figure  1.  The  reason 
no  distinction  between  contravarlant  and  covariant  vectors  occurs  for 
orthogonal  coordinates  can  be  seen  by  examining  (2-15).  When  orthogonal  unit 
base  vectors  are  used  the  1th  contravarlant  base  vector  defined  by  (2-15)  will 
be  equivalent  to  the  1th  covariant  base  vector. 


Contravarlant  and  covariant  vectors  are  special  cases  of  tensors.  The 
components  of  a  tensor  obey  the  following  transformation  rule: 
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(2-23) 

35 


The  exponent  N  of  the  Jacobian  Is  called  the  weight  of  the  tensor  field.  For 
N*0  the  tensor  field  Is  absolute;  If  N?*0  the  tensor  field  Is  relative  of 
weight  N.  The  tensor  is  contravarlant  of  order  r  and  covariant  of  order  s. 
By  setting  N*0  and  alternately  r»l,  8*1  the  previously  defined  transformation 
rules  for  contravarlant  (2-7)  and  covariant  (2-10)  vectors  are  recovered. 


II. 3.  The  Metric  Tensor 

In  mechanics,  the  notion  of  invariance  Is  of  fundamental  Importance.  Cer¬ 
tainly  the  meaning  of  a  physical  law  should  not  be  dependent  upon  the  choice 
of  the  coordinate  system;  therefore,  we  seek  to  formulate  the  equations  of 
fluid  mechanics  in  an  Invariant  form.  A  tensor  represents  an  Invariant 
quantity  whose  components  may  be  transformed  from  one  coordinate  system  to 
another  using  (2-23).  Thus,  equations  expressed  in  curvilinear  tensor  form 
are  Invariant  with  respect  to  coordinate  transformations. 


As  an  example  of  the  Idea  of  Invariance,  first  consider  a  scalar  function 
such  as  temperature.  A  scalar  la  an  absolute  tensor  of  zero  order.  The  tem¬ 
perature  at  a  point  on  a  sphere  has  the  same  value  no  matter  If  the  point  Is 
referenced  by  a  Cartesian  or  spherical  coordinate  system;  a  scalar  Is  Invar¬ 
iant.  By  the  same  token,  a  transformation  of  coordinates  only  changes  the 
components  of  a  vector;  the  vector  Itself  Is  an  Invariant.  A  vector  (contra- 
variant  components  and  covariant  basis) 


A 


(2-24) 


can  be  transformed  Into  an  alternate  coordinate  system  using  (2-7)  and  (2-10) 
to  obtain 


A  A^  9x^1  ^  rk-  „  jk-  _  7 

A  =  A  a  =  {A  -^Xa.  — “  A  a  6  =  A  a  »  A 


k 


(2-25) 


1  ax 

demonstrating  the  Invariance  of  the  vector. 

The  differential  arc  length  for  Cartesian  coordinates 


(ds)^  *  (dx)^  +  (dy)^  +  (dz)^ 


(2-26) 


Is  another  example  of  a  scalar  Invariant.  A  Euclidean  space  Is  one  In  which 
It  Is  possible  to  define  such  a  coordinate  system.  For  a  n-dlmenslonal 
Euclidean  space  the  differential  arc  length  Is 


(ds)' 


6  j^jdx^dx^ 


(2-27) 


where  the  Euclidean  metric  tensor  6 has  the  values 

lo  l^J 


(2-28) 


Now  the  x^  represent  Cartesian  coordinates.  Applying  a  coordinate  transforma¬ 
tion  to  dxl  and  dxJ  gives  the  general  elemental  arc  length 


3x^  ._k,^9x^ 


(ds)‘  -  6^,  dr)(^ 


'Ij 


dC*) 


3x 


(2-29a) 


where 


^Ik 


(2-29c) 


defines  the  metric  tensor.  Observe  that  the  metric  tensor  Is  symmetric, 
l.e.  =  gjj[.  Note  that  one  need  not  always  define  the  metric  tensor  as  a 

transformation  between  Cartesian  and  general  coordinates.  If  the  metric 
tensor  Is  known  In  one  coordinate  system  the  metric  tensor  In  any  other 
coordinate  system  can  be  found  using  the  tensor  transformation  rule  (2-23): 


’IJ 


laisnU 


I 

where  Is  the  metric  tensor  for  the  general  n  coordinates.  The  metric 
tensor  can  also  be  defined  by  considering  the  scalar  product  of  two  line 
elements  (2-13) 


(ds)‘ 


dr 


dr  ■  a,» 
~  ~1 


where 


®1J  -1* 


'j 


defines  the  metric  tensor.  Since  the  metric  tensor  Is  defined  In  the  same  way 
In  every  coordinate  system  It  is  an  Isotropic  tensor.  In  fact,  the  metric 
tensor  Is  the  most  general  second-order  Isotropic  tensor. 

An  Inverse  (contravarlant)  metric  tensor  g^^  Is  defined  by 


Jk 


(2-30) 


Notice  that  the  determinant  of  the  g^j  matrix  Is  the  square  of  the  Jacobian 
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For  two  dimensional  streamline-normal  coordinates,  (C  ,C  ) 
nonzero  components  of  the  metric  tensor  are 


(2-38) 


(s,n),  the 


.2  ^  n  ,2 

gii“  h 

o  (2-39) 

i.2  , 

822"  ^2  *  ^ 


where 

The 


is  the  radius  of  curvature  of  a  reference  streamline  or  body  surface, 
length  of  a  vector  Is  given  by 


|a|^  -  A.A  .  g^jA^A^ 


(2-40) 


from  which  the  associated  tensor  of  A^  is  defined  as 

Aj  -  gjjA^  (2-^0 

and  by  use  of  (2-30) 

A^  -  g^^A^  (2-42) 

This  provides  a  convenient  way  to  raise  or  lower  Indices.  Thus,  the  scalar 
product  of  two  vectors  can  be  expressed  as 

A  •  B  -  g..A^B^  -  A,B^  (2-43) 

~  ~  ®1J  J 

From  the  definition  of  a  scalar  product,  the  angle  between  two  vectors  Is 
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cos  6 


A  B 


and  the  angle  between  two  coordinate  lines  Is  given  by 


cos  6 


ij 


(2-44) 


(2-45) 


II. 4.  Physical  Coqionents  of  Tensors 

A  particularly  Inconvenient  property  of  curvilinear  tensors  Is  that  the 

dimensions  carried  by  the  tensor  components  do  not  necessarily  correspond  to 

the  physically  meaningful  dimensions  of  the  quantity  represented  by  the 

tensor.  For  example.  In  Cartesian  coordinates  each  coordinate  direction  has 

the  dimension  of  length.  In  contrast,  for  cylindrical  coordinates  not  every 

coordinate  direction  has  the  dimension  of  length.  Consider  the  contravarlant 

12  3 

velocity  components  In  cylindrical  coordinates  (^  ,5  )  =  (x,r,e); 


V 


1 


(2-46) 


It  Is  evident  that  does  not  possess  the  usual  (L/T]  dimensions  of  velocity; 
Instead,  Its  dimension  Is  [l/Tj.  Therefore,  tensor  components  may  not  repre¬ 
sent  physically  measurable  quantities. 

Methods  by  which  the  tensor  components  may  be  converted  to  physical  com¬ 
ponents  are  neatly  presented  by  Truesdell  (1953).  The  physical  component  can 
be  defined  as  the  magnitude  of  the  1th  component  projected  onto  the  1th  coor¬ 
dinate  direction: 


A(l)  “  (A^a. )  •  e,  (no  sum  on  1) 

~1  ~1 

■  /g  A^e  •  e 
*11  ~1  -1 

where  e^^  Is  a  unit  covariant  base  vector.  Based  on  this  definition  the 
physical  components  of  a  covariant  or  contravarlant  vector  are  given  by 
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A(l)  -  /gjj  A^  »  g^^  Aj  (no  sum  on  1)  (2-47) 

The  physical  components  of  higher-order  tensors  may  be  obtained  by  observing 
how  they  are  defined  In  terms  of  lower-order  tensors.  If 


c*  - 


then  an  equally  valid  equation  for  the  physical  components  Is 


(2-48) 


C(l)  -  A(1J)  B(J) 


(2-49) 


Using  (2-47)  In  (2-48)  yields 


cm.  1  B(jl 

^T~  ^  v'— 

**11 


(2-50) 


which  leads  to 


844  4 

A(lj)  -  /  Aj 

Further  examples  and  discussion  are  given  by  Truesdell. 


(2-51) 


II«5.  Derivatives  of  Tensors:  Covariant  Differentiation 

In  order  to  express  quantities  which  Involve  partial  derivatives,  for 
Instance  the  Laplaclan,  a  partial  derivative  operation  which  Is  valid  In  cur¬ 
vilinear  coordinates  Is  required.  The  need  for  a  general  partial  derivative 
Is  demonstrated  by  observing  that 


3A  -.1  .  3a. 

~  3A  .  .1  ~1 

— T  ■  — T  +  A  — r 
3?^  35^  3C^ 


(2-52) 


Immediately  It  Is  clear  that  the  variation  of  the  base  vectors  along  coordi¬ 
nate  lines  will  make  partial  differentiation  more  complicated  In  curvilinear 
coordinates.  In  Cartesian  coordinates,  the  partial  derivatives  always  give 
higher  order  tensors.  This  Is  not  true  for  curvilinear  coordinates.  For 


-14- 


m 


*•  - 


r/./. 


I  *  •  '  • 


.*•  i»V 


PR 


nri 


- 


1^ 


-  .  - 

'-y-'. .1 

.1 


where  l,j,k  are  In  cyclic  order. 

It  is  noted  that  both  J_  and  rj.  are  not  tensors;  however,  the 

covariant  derivative  (2-53)  does  preserve  the  tensor  transformation  rule  under 
coordinate  transformation  and  is  therefore  a  higher  order  tensor.  The  result 
of  equation  (2-53)  is  that  we  need  not  consider  explicitly  the  effect  of  coor¬ 
dinate  transformations  on  tensor  derivatives.  For  example,  since 

g  and  are  zero  In  Cartesian  coordinates,  they  are  zero  in  all  coordl- 
*ij,k  *,k 

nate  systems. 

From  equation  (2-53),  the  covariant  derivative  of  a  covariant  vector  A 
with  respect  to  is  given  by 


“"l  k 


(2-57) 


Similarly,  the  covariant  derivative  of  a  contravariant  vector  A*  is  given  by 


.1  .  a'" 

3,J 


(2-58) 


The  covariant  derivative  of  the  physical  components  of  a  contravariant  vector 
may  be  found  by  introducing  the  physical  components  into  (2-58) 


where 


A^  «  g*^^  -t  D(im)  A(l)  +  r  (Imn)  A(n)}  (2-59) 

*,m  ®11  *mm  v  /  '  ^ 


3  .  „-l/2  _L. 

3C(m)  *mm  ^^m 


(2-59a) 
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(2-59b) 


and  r(lmn)  is  the  physical  component  of  the  Chrlstoffel  symbol  defined  by 


^  1/2  -1/2  -1/2  1 
r(imn)  -  gjj 


(2-59c) 


If  Al  «  a1  (£(t),t),  then  the  intrinsic  (total)  derivative  of  A*  with  respect 


to  time  is 


iir  *  ^ 
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Since 

at 

of  A  . 


is  the  velocity. 


this  Is  an  expression  for  the  material  derivative 


With  the  foregoing  tensor  concepts,  we  are  now  In  a  position  to  recast 
many  of  the  familiar  operations  of  vector  analysis  Into  curvilinear  tensor 
form.  Although  these  are  not  directly  employed  In  subsequent  sections,  they 
are  summarized  In  the  Appendix  since  they  are  useful  In  gaining  Insights  Into 
the  equations  of  motion  and  relating  the  tensor  equations  to  the  more  familiar 
vector  forms. 


UATIOMS  OF  FLUID  MECHANICS 


The  equations  governing  fluid  motion  In  a  nonorthogonal  coordinate  system 
can  be  derived  using  the  methods  of  curvilinear  tensor  analysis.  The  results 
of  the  previous  section  provide  the  necessary  framework  for  writing  these 
equations  In  an  Invariant  form;  the  physical  meaning  of  the  equations  will  be 
consistent  under  any  coordinate  transformation.  Alternatively,  the  equations 
could  be  derived  by  the  methods  of  vector  analysis.  The  details  of  the  vector 
analysis  approach  are  presented  by  Stern  and  Yoo  (1986). 

The  majority  of  equations  are  given  using  tensor  components 
(contravarlant  or  covariant)  which  do  not  necessarily  carry  the  proper 
physical  dimensions  of  the  quantity  they  represent.  In  a  few  example  cases 
the  equations  are  written  substituting  the  physical  components  for  the  tensor 
components.  As  a  rule,  tensor  components  are  more  easily  manipulated  and  we 
prefer  their  use  over  physical  components  when  In  the  process  of  deriving 
general  results. 

Many  of  the  equations  given  here  are  straightforward  generalizations  of 
their  Cartesian  tensor  counterparts.  The  reader  Is  encouraged  to  compare  the 
equations  presented  here  to  those  In  Cartesian  tensor  form  (see,  for  example, 
Currie,  1974).  Loosely  speaking,  the  Cartesian  tensor  form  of  the  equations 
can  be  generalized  by  replacing  partial  derivatives  with  covariant  deriva¬ 
tives,  making  the  velocity  a  contravarlant  vector,  and  keeping  the  proper 
overall  covariant  and/or  contravarlant  order. 

The  material  derivative  represents  the  rate  of  change  of  a  quantity  as 
seen  by  an  observer  moving  with  a  fluid  particle.  Therefore,  the  material 
derivative  of  a  general  tensor,  which  may  be  a  function  of  both  space  and 
time,  is  given  by  Its  Intrinsic  (total)  derivative 


M  lA  .  A  ^ 
Dt  “  3t  ^,J  dt 


(3-1) 


Now  represents  the  jth  contravarlant  velocity  component.  Letting  V-^ 


the  expression  for  the  material  derivative  becomes 


(3-2) 


M  -  M  +  A  vJ 

Dt  3t 

Notice  that  this  expression  Is  quite  similar  to  the  material  derivative 
defined  using  Cartesian  tensors,  but  now  the  usual  partial  derivative  has  been 
replaced  by  the  covariant  derivative.  If  the  general  tensor  A  represents  a 
velocity,  equation  (3-1)  may  be  written  In  terms  of  contravarlant  velocity 
components  as 


If  . 


(3-3a) 


(3-3b) 


This  represents  the  fluid  acceleration  in  general  tensor  form. 

The  material  derivative  of  a  volume  Integral  Is  given  by  the  Reynolds 
Transport  Theorem.  The  theorem  Is  useful  In  deriving  Eulerlan  conservation 
laws  from  their  Lagranglan  equivalents.  Generalizing  the  Cartesian  tensor 
form  of  the  theorem  yields 


f  A  dV  «  /  {-ff  AV,J}  dV 


(3-A) 


where  A  Is  a  general  tensor  and  V  represents  a  material  volume.  An  alternate 
and  perhaps  more  familiar  arrangement  Is 


f  /  AdV=/  {|f  +  (AV^)^^}  dV 


(3-5) 


III.l.  Equation  of  Continuity 

The  conservation  of  mass  or  continuity  equation  may  be  derived  by  letting 
the  general  tensor  A  represent  the  fluid  density  p  In  (3-4),  by  which 


(3-6) 


Since  the  volume  considered  Is  arbitrary  the  integrand  must  be  equal  to  zero: 
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(3-7) 


If  the  fluid  Is  Incompressible  the  volume  of  a  material  element  must  remain 
constant.  In  other  words 


ll-o 


(3-8) 


Hence,  the  statement  of  mass  conservation  for  an  Incompressible  fluid  becomes 


v;..o 


(3-9) 


which  Is  V  •  V  ■  0  In  vector  notation. 


A  useful  corollary  to  the  Reynolds  Transport  Theorem,  derived  using  the 
conservation  of  mass  statement.  Is 


It  ^  PA  dV  >  /  P  ^  dv 


(3-10) 


Substitution  of  the  physical  velocity  components  for  the  contravarlant 
components  In  (3-9)  gives 

a  -1/2 

-1/2  3V(1)  _  ®li  .  _1  -1/2 .  ,,  ..V 

«11  77l  V(i)  — r  +  «kk  V(k)  *  0  (3-11) 


Note  the  sum  over  both  1  and  k.  Using  (2-55b) 
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and  the  continuity  equation  may  be  written  In  the  form 

1  _i_  .  0 

72  ,,1  ^~T72  ^ 


(3-12) 


’ll 


III. 2.  Hoaentum  (Navler-Stokes)  Equations 

The  conservation  of  linear  momentum  (Newton's  Second  Law)  for  an 
arbitrary  material  volume  can  be  expressed  as 
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— /  p  V*AdV-/  pF  •  i  d¥  +  f  t  •  Z  dS  (3-13) 
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where  V,  F,  and  t  are,  respectively,  the  velocity,  body  force,  and  surface 
stress  vectors.  Because  the  absolute  direction  of  the  base  vectors  Is,  In 
general,  variable,  each  force  should  be  resolved  along  a  constant  direction 
Z.  The  need  for  this  refinement  can  be  seen  by  considering  the  forces  on  a 
fluid  element  defined  by  a  cylindrical  coordinate  system. 

The  nature  of  the  surface  forces  can  be  deduced  by  examining  the  forces 
acting  on  an  elemental  volume  of  fluid.  If  the  volume  in  (3-13)  is  decreased 
to  a  point  while  preserving  its  shape,  we  find  that  the  surface  forces  are  in 
a  state  of  local  equilibrium: 

11m  t  •  £  dS  =  0  (3-14) 

¥*■0  r  ^ 

where  r  is  a  characteristic  dimension  of  the  volume.  We  need  to  find  how  the 
stress  on  an  arbitrary  surface  is  related  to  the  stresses  on  the  coordinate 
surfaces.  Figure  1  shows  a  typical  volume  element  and  the  stress  vectors 
acting  on  each  surface.  Since  the  forces  are  in  local  equilibrium  we  have 


t  dS  »  t,  dS,  +  t_  dS-  +  t_  dS_ 

~n  ~1  1  ~2  2  ~3  3 

where  dS  refers  to  the  area  of  the  arbitrary  surface  with  normal  n. 
the  normal  vector  is  composed  of 

1  1 
n  “  n  a,  “  n,  a 
~  ~1  i  ~ 


(3-15) 


Note  that 


(3-16) 


The  outward  normal  vector  (not  necessarily  a  unit  vector)  for  each  coordinate 
surface  is  a^.  The  surface  area  of  each  face  of  the  element  is  related  to  dS 
by 


dS,  =  a,  •  n  dS 
i  ~1  ~ 


Si 


(Oj  8“^)  dS  =  n^  dS 


(3-17) 


Substituting  for  each  dSj  in  (3-15)  yields 


The  stress  vector  for  each  coordinate  surface  Is  then 


(3-18) 


where  o 

Now  t 
~n 


defines  the  stress  tensor  which  represents  the  1th  component  of  t 
can  be  written  as 
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~n 


(3-19) 


Evidently,  the  1th  contravarlant  component  of  stress  on  the  surface  with 
outward  normal  n  Is 


t^  -  o^^nj  (3-20) 

The  stress  tensor  components  are  shown  In  figure  2.  Observe  that  does  not 
necessarily  represent  a  normal  stress;  it  is  simply  the  stress  In  the  1th 
coordinate  direction.  Only  In  an  orthogonal  coordinate  system  will  be  a 
normal  stress.  Since  we  are  considering  fluids  which  cannot  have  Internal 
stress  couples  (moments),  angular  momentum  can  be  conserved  only  if  the  stress 
tensor  Is  symmetric,  i.e. 


IJ  Ji 

This  result  can  be  derived  by  considering  the  moment  of  momentum  equation. 

Utilizing  the  expression  for  the  stress  on  an  arbitrary  surface  (3-20), 
the  momentum  equation  may  be  rewritten  as 

/y  P  ^  ll  dV  -  dV  (3-21) 

where  (3-10)  and  (A-16)  have  been  used.  Since  the  direction  t  is 
constant  (£  *  0)  and  both  the  volume  and  1  are  completely  arbitrary,  the 

1  ,  K 

equality  can  be  satisfied  only  If 


This  result  Is  known  as  Cauchy's  equation  of  motion  and  It  Is  valid  for  any 
continuum.  In  order  to  solve  this  equation  a  relation  between  the  stress 
tensor  and  the  velocity  field  Is  required. 

Although  a  complete  development  of  the  constitutive  equation  for  a 
Newtonian  fluid  Is  beyond  the  scope  of  this  report,  a  few  details  of  the 
derivation  will  be  shown  In  order  to  make  the  origin  and  form  of  the  final 
result  clear.  A  rigorous  derivation,  using  Cartesian  tensors.  Is  available  In 
many  texts  (e.g.  Batchelor  (1967)). 

In  a  motionless  fluid  the  force  on  a  fluid  element  Is  due  solely  to  the 
Internal  pressure.  For  an  Incompressible  fluid,  this  Is  the  mechanical  pres¬ 
sure.  If  the  fluid  Is  compressible,  the  pressure  would  be  the  thermodynamic 
pressure  specified  by  an  equation  of  state.  For  a  static  fluid  the  stress 
tensor  Is  Isotropic  and  can  be  written  as 

*  —  p  (3—23) 

From  this  the  stress  on  an  elemental  area  can  be  determined  using  (3-19): 

/  ^  /lx 

“  (-Pg  nj)  a^  -  (-pn  )  a^  =  -  pn 

Consequently,  the  pressure  acts  normal  to  the  surface. 

For  a  moving  fluid,  a  more  general  form  of  the  stress  tensor  should  be 
considered: 


.  -  pgiJ  +  (3-24) 

The  additional  stress  due  solely  to  the  motion  of  the  fluid  Is  the  devlatoric 
stress  The  devlatoric  stress  Is  assumed  to  be  proportional  to  the  fluid 


deformation  rate 


and  V 


k 

.J 


can  be  decomposed  Into  Its  symmetric  and  antisymmetric  parts 


v'^1  =  (eij  +  fiij)  (3-25) 

where  e^j  =  ~  (^l,j  ^IJ  “  T  1^*  Since  the  stress 
tensor  must  be  symmetric  the  devlatorlc  stress  cannot  depend  upon  the  antisym¬ 
metric  portion  of  the  deformation  rate.  Assuming  that  a  linear  relation 
between  the  devlatorlc  stress  and  strain  rate  (eij)  exists. 


T 


IJ 


(3-26) 


1  4J^  gj 

where  D  Is  a  fourth-order  Isotropic  tensor  (In  fact.  It  Is  a  composition 
J^QI  1.  i 

of  g  ),  Hence,  t  Is  related  to  e,  by  a  linear  combination  of  Isotropic, 

tm 

fourth-order  tensors: 


Ij  ,  ij  £m  ,  .  ml  tj .  11  jm. 

^  «  ®tm  ^  ®lm 


(3-27) 


where  u  Is  the  dynamic  viscosity  and  X  Is  the  second  viscosity  coefficient. 
This  can  be  simplified  to 


.  Xg^J  +  2P  e^J 


(3-28) 


where  =  g^*^  g^™  ®£in*  Incompressible  fluid  (V^^-O)  the 

Newtonian  stress  tensor  Is 


a  -  -  -  g  “p  +  2u  e^^  -  -g^^P  +  '^1  ,k^  (3-29) 


IJ 


IJ. 


In  (3-22)  the  divergence  of  the  stress  tensor  Is  required.  Carrying  out  this 
operation,  using  the  fact  that  g^^  »•  0,  gives 


a 


iJ 

,J 


-g‘Jp 


.i 


(3-30) 


For  most  fluid  mechanics  applications  It  can  be  assumed  that  the  space  Is 

Euclidean.  This  guarantees  that  yJ  .  yJ  , which,  by  the  equation  of 

,  k  J  ,  J  k 

continuity.  Is  equal  to  zero.  Then  the  divergence  of  the  stress  tensor  Is 
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Notice  that  the  only  free  Index  Is  1,  which  makes  the  result  a  contravarlant 
vector. 

The  Navler-Stokes  equations  for  an  Incompressible  fluid  are  obtained  by 
substituting  (3-31)  Into  the  equation  of  motion  (3-22).  The  result  Is 


If  * 


L™  +  iL  e'xnvl  +  pi 

’  p  p  ®  ,mn 


(3-32) 


Replacing  by  V(l)  and  making  use  of  (2-59)  and  (A-15)  yields  the 

following  Navler-Stokes  equations  with  physical  components  as  dependent 
variables: 


+  V(m)  +  D(lm)  V(m)  V(l)  +  r(lmn)  V(m)  V(n)  +  g(lm)  3j(iO 

-v{V^V(l)  +  2E(lm)  2C(imn)  +  P(l)  V(l)  +  Q(ln)  V(n)}  -  F(l)  =  0 


where 


F(l)  - 


(3-33) 

(3-3A) 


and  E(lm),  C(lmn),  P(l)  and  Q(ln)  are  defined  In  equations  (A-15a)  through  (A- 
15d).  The  Laplaclan  of  V(i)  is  given  by  equation  (A-14)  as 


3C  35 


(3-35) 


These  equations  are  rather  cumbersome  when  expanded  for  a  given  value  of  1 


because  many  terms  include  Implied  summations.  Even  though  both  g(lj)  and 
r(ljk)  are  symmetric  they  have  6  and  18  Independent  components,  respectively. 
Note  that  using  the  physical  velocity  components  has  Introduced  additional 
terms  like  D(lm)V(m)V(l)  Into  the  equation  which  are  not  present  In  the 
contravarlant  form  (3-32). 


III. 3.  Energy  Equation 


By  applying  the  first  law  of  thermodynamics  to  a  fluid  volume,  an 
equation  expressing  the  conservation  of  energy  can  be  obtained.  Since  energy 
Is  a  scalar,  the  formulation  of  the  energy  equation  provides  a  general  example 
of  how  transport  equations  for  other  scalar  quantities  may  be  derived.  The 
first  law  states  that  for  a  given  system  the  rate  of  change  of  total  energy  Is 
equal  to  the  rate  of  work  done  plus  the  rate  of  heat  added.  For  an  arbitrary 
fluid  volume  we  have 

P(e+^y*V)dV-  /  y*tdS+/  V*pFdV-  /  £  •  n  dS  (3-36) 

where  e  Is  the  Internal  energy  per  unit  mass  and  £  denotes  the  conductive  heat 
transfer  (positive  for  heat  leaving  the  volume).  The  left  hand  side 
represents  the  rate  of  change  of  total  (Internal  and  kinetic)  energy  and  the 
right  hand  side  Is  composed  of  the  rate  of  work  done  by  surface  forces,  rate 
of  work  done  by  body  forces,  and  the  rate  of  heat  addition  respectively.  Con¬ 
verting  (3-36)  to  general  tensor  form  and  using  (3-10)  yields 

/  0  ^  j  g,,vV)  d*  -  /  s  tV  dS  +  /  eg  v^p'd*  -  /  g  dV  dS 

*  S  J  ^  S  J 

After  substituting  (3-20)  for  the  surface  force  vector  and  using  the 

divergence  theorem  on  the  surface  Integrals,  we  obtain 

/  P  (e  +  i  dV  «  /  ((a  +  p  -  q*^  ]dV  (3-38) 

Since  the  volume  considered  is  completely  arbitrary  the  equality  holds  only  if 

^  T  P'^1^^  "  ‘’1^1  (3-39) 

The  kinetic  energy  term  may  be  rewritten  as  follows: 
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—  vV) 

Dt  '^2 


1  ,1  ^ 

2  Dt 


g 

*1J  Dt 


j 


DV. 
1  Dt 


J 


+  I  g 


IJ  ^  Dt 


The  mechanical  energy  equation  (Cauchy's  equation  (3-22)  multiplied  by  V^)  may 
now  be  Identified  and  subtracted  out  from  (3-39),  leaving 


P 


Dt 


(3-40) 


which  Is  a  transport  equation  for  the  Internal  energy.  The  heat  flux  vector 
can  be  related  to  the  temperature  field  through  Fourier's  law: 


q^  -  -  g^^KT  ,  (3-41) 

where  K  Is  the  thermal  conductivity.  Substituting  for  the  heat  flux  in  (3-40) 
gives 

P  §1-  (g‘JKT,,)_,  (3-42) 

The  last  term  of  the  above  equation  represents  the  increase  In  Internal 
energy  due  to  the  work  done  In  deforming  a  fluid  element;  In  other  words, 
conversion  of  mechanical  energy  to  thermal  energy  by  the  action  of  surface 
forces.  Replacing  the  stress  tensor  with  Its  incompressible  fluid  form  (3- 
29),  we  have 


-  (-pg^^  +  2pe^^) 

-  -pV^,  +  2Me^\.  . 

»J  ^  >  J 


-  2ue^^  V 


IJ 


(3-43) 


Further  Insight  can  be  gained  by  splitting  Vj  j  Into  Its  symmetric  and  anti¬ 
symmetric  parts: 


®lj  ■*■  "^IJ 


fSS^ 


Now  substituting  this  result  Into  (3-43)  gives 

.  2pe^J  (e^^ 

and  recognizing  that  the  multiplication  of  a  symmetric  tensor  and  anti¬ 
symmetric  Censor  Is  zero,  we  obtain 


4  ■  j  =  2pe^^ej^j 


(3-44) 


•  w  * 


where  4  denotes  Che  dissipation  function.  The  dissipation  function  Is  a 
positive  definite  quantity  which  represents  the  Irreversible  conversion  of 
mechanical  energy  Into  thermal  energy  due  to  the  action  of  viscosity.  The 
energy  equation  for  an  Incompressible  fluid  may  now  be  written  as 


<f  kt  ,)  ,  +  » 


(3-45) 


III. 4.  CaaCloQ  on 


jroxlmatioas 


At  this  point  It  seems  appropriate  to  make  some  general  comments 

regarding  the  Invariance  of  the  governing  equations.  Recall  that  a  tensor 

represents  a  quantity  which  Is  Invariant  under  coordinate  transformations. 

The  equations  presented  above  are  written  In  general  tensor  form  and 

therefore,  are  Invariant.  However,  when  the  covariant  derivatives  are 

expanded  and  physical  components  are  Introduced,  some  parts  of  the  equations 

1  k 

will  not  be  tensors.  For  example,  f..  A  Is  a  part  of  the  covariant 

j  k 

derivative,  but  It  Is  not  a  tensor.  The  covariant  derivative  Itself  Is  a 
censor  and  thus  represents  an  Invariant  quantity.  It  Is  also  Important  to 
note  that  after  expansion  some  Individual  terms  of  an  equation  may  appear  to 
be  Identical,  but  In  reality  they  may  have  originated  from  different  Invariant 
groups  Chat  do  not  have  the  same  physical  meaning. 


-28- 


Another  Important  point  to  be  noted  Is  that  when  approximations  are  made 
to  an  Initially  Invariant  equation  the  resulting  equation  may  lose  the  invar¬ 
iance  property.  A  familiar  example  of  this  Is  provided  by  the  approximation 
of  the  Navler-Stokes  equations  to  form  the  boundary-layer  equations.  Since 
the  boundary-layer  approximation  neglects  only  certain  terms  in  the  original 
invariant  groups,  the  resulting  equations  are  not  Independent  of  the 
coordinate  system.  The  implications  of  this  fact  were  discussed  some  time  ago 
by  Kaplun  (1954). 

Another  problem  Inherent  In  making  approximations  based  on  the  order  of 
magnitude  of  Individual  terms  Is  that  of  estimating  the  contribution  from  the 
geometric  coefficients.  Clearly,  one  should  consider  for  each  terra  the  order 
of  magnitude  of  both  the  Independent  and  dependent  variables.  Even  though  a 
quantity  composed  of  dependent  variables  may  be  of  small  order.  It  may  in  turn 
be  multiplied  by  a  large  geometric  coefficient,  such  that  the  net  result  Is  a 
terra  of  first  order  significance. 


IV.  EQDATIOMS  FOR  TDRJMTLEHT  FLOW 


It  Is  generally  assumed  that  the  e(;uations  derived  In  the  previous  sec¬ 
tion  apply  in  a  turbulent  flow  at  all  instants  of  time.  For  most  practical 
purposes,  however,  we  follow  Reynolds  and  decompose  the  instantaneous  values 
of  flow  variables,  such  as  velocity  and  pressure,  into  a  mean  value  and  a 
fluctuation  about  that  mean.  Henceforth,  the  Instantaneous  velocity  is 

J  J 

decomposed  into  +  v  and  the  Instantaneous  pressure  is  decomposed  into  p  + 
p'.  The  equations  governing  the  mean  flow  are  then  derived  from  those  given 
in  the  previous  section  by  making  this  substitution  and  taking  the  mean 
values. 

IV. 1.  Continuity  Equation 

Introducing  the  mean  and  fluctuating  components  into  the  continuity 
equation  (3-9)  gives 

V^.  +  v^  =  0  (4-1) 

and  averaging  the  equation  we  have 


> 


1 


*  0 


From  (3-12),  this  may  be  written  in  conservative  form  as 


sc 


V(i)}  »  0 


or,  in  non-conservative  form,  as 


3V(i) 

^tnr 


+  iD(li)  +  r(nnl)J  V(l) 


0 


(4-2) 


(4-3) 


(4-4) 


Thus,  the  mean  velocity  field  satisfies  the  usual  continuity  equation  (3-9)  or 
(3-12).  The  fluctuating  velocity  field  also  satisfies  the  continuity  equation 
since  subtracting  (4-2)  from  (4-1)  yields 
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IV.2.  MDaentua  (Reynolds)  Equations 


Substituting  the  mean  and  fluctuating  components  of  velocity  and  pressure 
Into  the  Navler-Stokes  equations  (3-32)  and  Reynolds  averaging  the  result 
gives 


If  * 


-  —  g^™p  +  vg"*"  -  (v^v“) 

p  *  ,mn  ,m 


+  F 


(4-6) 


These  are  the  momentum  equations,  or  Reynolds  equations,  for  turbulent  flow. 
These  equations  are  Identical  to  the  original  Navler-Stokes  equations  except 
for  the  appearance  of  the  Reynolds  stress  tensor  v^v*”.  The  additional  Reynolds 
stress  terms  arise  as  a  result  of  the  nonlinearity  of  the  equation  of  motion, 
and  can  be  expressed  In  terms  of  their  physical  components  as 

(v^v")  =>  tD(lm)  +  D(mm)  +  r(nnm)]  v(l)v(m) 

,ni  11  oC  (.ni) 

+  r(lmn)  v(m)v(n)}  (no  sum  on  1)  (4-7) 


Therefore,  by  utilizing  equations  (3-33),  the  Reynolds  equations  for  turbulent 
flow  become 


3V(1) 

3t 


V(m)  +  D(lm)  V(m)  V(l)  +  r(lmn)  V(m)  V(n)  + 


[D(lm)  +  D(mm)  +  r(nnm)]  v(l)  v(m)  +  r(lmn)  v(m)  v(n)  = 


■  ''^^^V(l)  +  2E(lm)  +  2C(lmn) 


3V(n) 

3C(m) 


+  P(l)  V(l)  +  Q(ln)  V(n)}  -  F(l)  =  0  (no  sum  on  1)  (4-8) 


Due  to  the  presence  of  the  six  unknown  Reynolds  stresses  v(l)v(m) ,  the  contin¬ 
uity  and  three  momentum  equations  (4-8)  are  no  longer  sufficient  to  determine 
the  mean  flow.  It  Is,  therefore,  necessary  to  furnish  expressions  for  the 
unknown  Reynolds  stresses  so  that  a  closure  of  the  equations  can  be  achieved. 


IV.3.  Equations  for  the  Reynolds  Stresses 

In  addition  to  the  Reynolds  equations  (4-8)  for  the  mean  motion,  a  paral¬ 
lel  set  of  equations  which  describe  the  transport  of  Reynolds  stresses  can  be 
derived  using  the  following  scheme: 


v^M^ 


where  represents  the  1th  Instantaneous  momentum  equation.  The  resulting 
transport  equations  for  the  Reynolds  stresses  are 


3  /  f  j\  .  .,™/  i  .  j  ®  ..1  f  ™  i  j  ™\ 

^  (v  V^)  +  V  (v  V  )  „  +  V  V  V  +  V  V  V'’^  +  (v  V  V  ) 
dt  ,m  ,m  ,m 


,  jm  ,v^p\  Im  ,v^p*.,  p'  .  Im  1  jm  1  . 

=  -  g  (“T^)  +  g  „  (g  +  g-^  V  ) 

P  °  p  .m  p  .m  .m 


,  1  i  j  1  V 

+  ug  (v  V  +  V  V  ) 
,mn  ,mn 


(4-9) 


The  first  two  terms  make  up  the  material  derivative  (time  and  convective  deri¬ 
vatives)  following  the  mean  flow.  The  following  two  terms  represent  the 
production  of  v^v^  by  the  rate  of  strain  of  the  fluid.  The  covariant  deriva¬ 
tive  of  the  triple  correlation  terms  stands  for  turbulent  diffusion  resulting 
from  the  convective  motion  of  the  turbulence.  The  first  term  on  the  right 
hand  side  of  the  equation  is  turbulent  diffusion  due  to  pressure  fluctuations. 
It  Is  often  grouped  together  with  the  triple  velocity  correlation  term.  The 
next  term  Is  usually  called  the  "pressure-strain"  or  "tendency-towards-lso- 
tropy"  term.  It  leads  to  the  redistribution  of  energy  among  the  six  Reynolds 
stresses  v^v^  and  an  approach  towards  the  statistically  most  probable  state. 
In  which  the  turbulence  Is  Isotropic,  The  last  term  Involving  the  viscosity 
represents  the  transport  of  energy  by  molecular  action  (viscous  diffusion)  as 
well  as  Its  dissipation  Into  heat. 

If  we  expand  the  covariant  derivatives  and  Introduce  the  physical  compo¬ 
nents  of  velocity  and  geometric  coefficients,  the  transport  equations  (4-9) 
can  be  expressed  as 
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+  [D(lm)  +  D(jm)I  V(m) 


ft  [v(l)v(j)]  +  V(m)  [v(l)v(j)j 


+  r(lmn)  V(m)  v(j)v(n)  +  r(jmn)V(m)  v(l)v(n) 


+  v(j)v(m)  [||^  +  D(lm)  V(l)  +  r(lmn)  V(n)  ] 

+  v(i)v(m)  +  D(jra)  V(j)  +  r(Jmn)  V(n)] 

+  [v(l)v(j )v(m)]  +  [D(lm)  +  D(jm)  +  D(mm)  +  r(nnm)]  v(l)v( j )v(m) 


+  r(lmn)  v(j)v(m)v(n)  +  r(jnin)  v(l)v(ni)v(n) 


-  lE(lJ)  ECU)  ^^1  -  IC(lJn)  *  CUin)) 


*  ^  Its  *  US'  *  ^  -  «J‘> 


+  [C(ljn)  t  C<llii)I 


+  v{.<I)  V^v(j)  +  V<J)  V^v(l) 


t  2E(Jm)  v<l)  2E(lm)  .(J) 


+  2C<lmti)  v(J)  t  2C(jinn)  v(l) 

+  (P(i)  +  P(j)]  v(l)v(j )  +  Q(ln)  vdj )v(n)  +  Q(jn)  v(l)v(n)  }  (no  sum  on  1  or  j ) 

(4-10) 


The  physical  Interpretation  of  the  various  terms  in  the  above  equation  Is  much 
the  same  as  It  was  for  equation  (4-9),  For  clarity,  however,  we  shall  examine 
some  of  the  terms  in  greater  detail.  It  is  clear  that  the  first  term  Is  the 
time  derivative  of  v(l)v(j),  and  the  next  four  terms  represent  the  convective 
derivatives  which  describe  the  transport  of  v(l)v(j)  by  the  mean  flow.  The 
following  two  terms  again  represent  the  production  of  v^l)v(j )  by  the  mean 


race  of  strain  In  the  fluid.  It  Is  noted  that  some  quantities  In  the  produc¬ 
tion  terms  have  exactly  the  same  form  as  those  which  appeared  In  the 
:onvectlve  terms  although  their  origin  Is  quite  different.  In  the  modeling  of 
Reynolds  stress  equations,  It  Is  Important  to  use  the  correct  expression  for 
the  production  term  In  order  to  retain  the  physical  significance  of  the 
modeled  terms.  As  mentioned  before,  terms  Involving  triple  velocity  cor¬ 
relations  represent  turbulent  diffusion  due  to  the  velocity  fluctuations. 

On  the  right  hand  side  of  equation  (4-10),  the  first  three  terms  make  up 
turbulent  diffusion  due  to  the  pressure  fluctuations,  while  the  following 
three  terms  Involving  pressure  fluctuations  are  the  pressure-strain  terms.  It 
Is  Important  to  note  that  all  pressure-velocity  correlations  without  deriva¬ 
tives  will  cancel  upon  summation.  However,  the  first  term  alone  does  not  rep¬ 
resent  the  total  turbulent  diffusion  due  to  the  pressure  fluctuations. 
Similarly,  Che  fourth  term  Itself  does  not  make  up  Che  overall  pressure-strain 
term.  The  pressure-strain  terms  In  general  curvilinear  coordinates  do  not 
possess  a  clear  physical  significance  similar  to  that  described  In  Nash  and 
Patel  (1972)  for  orthogonal  coordinates  since  v(l)v(l)  (unsummed)  no  longer 
represents  the  turbulent  stresses  normal  to  the  surfaces  of  a  fluid  element. 
The  redistribution  of  energy  and  approach  Cowards  Isotropy  should  be 
Interpreted  by  examining  the  turbulent  kinetic  energy,  which  Is  a  contraction 
of  the  Reynolds  stresses. 

It  Is  seen  that  the  transport  equations  (4-10)  contain  many  additional 
unknown  quantities  such  as  triple  correlations  and  pressure-strain  terms. 
Approximations  must  be  made  either  at  this  or  ocher  levels  for  the  unknown 
correlation  terms  In  order  to  obtain  a  closure  of  the  governing  transport 
equations. 


IV.4.  ItarbulenC  Kinetic  Energy  Equation 

A  quantity  often  used  In  turbulence  closure  schemes  Is  the  turbulent 
kinetic  energy  defined  by 


or 


k 


1 

Y  V*  V 


11  1  i  j 

2  ''  '"l  “  2  '^IJ 


k  »  Y  cos  v(l)v(J) 


(4-11) 
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(4-12) 


A  transport  equation  for  turbulent  kinetic  energy  can  be  derived  by  contrac¬ 
ting  (multiplying  by  -i-  gj^ j  )  equation  (4-9): 

3  k  m,  j  m  1  ,1  1  j  mn 

yt  +  V  +  g^j  V  V  +  (y  J?ij  VV  V  +  (-^)^„  -  vg  k^^^ 


mn  1  1  - 

+  vg4  4g  V  v-’^  =  0 


(4-13) 


The  first  two  terms  represent  the  material  derivative  following  the  mean  mo¬ 
tion.  The  next  term  stands  for  the  production  of  turbulent  kinetic  energy  by 
the  mean  rate  of  strain,  or  the  extraction  of  mean  flow  energy  by  the  tur¬ 
bulence.  The  next  three  terms  represent  diffusion  due  to  turbulent  velocity 
fluctuations,  pressure  fluctuations  and  molecular  action,  respectively.  The 
last  term  Involving  viscosity  represents  the  decrease  of  turbulent  kinetic 
energy  due  to  viscous  dissipation.  Note  that  this  Is  equal  to  the  true  dis¬ 
sipation  (3-44)  only  for  homogeneous  turbulence. 

Note  that  the  contraction  of  the  pressure-strain  term  Is  Identical  to 
zero,  by  virtue  of  the  equation  of  continuity,  even  In  a  general  curvilinear 
coordinate  system.  In  other  words,  there  Is  no  net  contribution  from  the 
pressure-strain  terms  to  the  turbulent  kinetic  energy.  This  clearly  suggests 
that  the  main  role  of  the  pressure-strain  terra  Is  In  the  redistribution  of 
energy  towards  Isotropy. 

When  physical  components  of  and  v^  are  used,  the  transport  equation 
for  k  becomes 

It  affmT  ®lj  v(J)v(m)  V(l,m)  +  v(m)  Ij  cos  0jj  v(l)v(j)] 


3C(m)  1  +  iD(mm)  +  r(nnm)]  - vv\  +  e  -  0 


(4-14) 


where  the  last  term  Is  the  viscous  dissipation 


€  »  V  cos  0^jg(mn)  v(l ,m)v(.1  ,n] 


(4-15) 


v'V 


and  V(l,m)  and  v(l,m)  are  the  physical  components  of  the  covariant  derivative 

and  v^  ,  i.e., 

,m  ,m*  ' 


V(i.m)  -  yi^  .  |vm  ^  ^  j.(ims)  VCs) 


(4-16) 


Although  it  is  possible  to  group  the  turbulent  diffusion  terms  due  to  the 
pressure  fluctuations  in  the  following  form  using  the  equation  of  continuity 
(4-4) : 


3  ,v(m)p',  p’  3v(m) 


(4-17) 
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the  second  term,  however,  does  not  represent  the  pressure-strain  term.  In  the 
turbulence  modeling  of  the  k  equation,  it  is  therefore  Incorrect  to  model 
these  two  terms  separately.  Similarly,  the  triple  velocity  correlation  term 
in  equation  (4-14)  should  not  be  split  into 


3v(m) 


3^^  [j  cos  v(l)  v(J)  v(m)]  -  j  cos  6^^  v(l)  v(J) 


(4-18) 


and  modeled  Independently  since  both  terms  represent  only  part  of  the  turbu¬ 
lent  diffusion  due  to  the  velocity  fluctuations. 

IV. 5.  Equation  for  Viscous  Dissipation  e 

A  transport  equation  for  the  dissipation  term  e  defined  in  (4-15)  can  be 
derived  using  the  following  scheme: 


'^^Ij  «!r  ^!r 


(4-19) 


where,  as  before,  represents  the  1th  instantaneous  momentum  equation. 


The  resulting  transport  equation  is 
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It  Is  noted  again  that  the  transport  equation  for  e  Involves  many  additional 
unknovm  quantities  which  can  not  be  expressed  In  terms  of  lower-order  correla¬ 
tion  by  simply  deriving  their  corresponding  transport  equations.  In  the  well 
known  k-e  turbulence  model,  all  terms  on  the  right  hand  side  of  equation  (4- 
20),  except  the  molecular  diffusion  term  vg  ^  nin  ”  modeled  to 
obtain  a  closure  equation. 


The  transport  equation  for  e  can  be  written  In  terms  of  the  physical 
components  by  utilizing  the  expressions  given  In  equation  (2-53)  for  covariant 
derivatives  and  the  usual  definitions  for  physical  components  of  tensor  quan¬ 
tities.  No  such  attempt  has  been  made  here,  not  only  due  to  the  obvious 
complexity  of  the  result  but  also  because  most  closure  schemes  use  only 
modeled  forms  of  this  equation. 
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V.  EDDY-VISCOSITY 


k-e  HOOBL 


The  classical  Bousslnesq  eddy-vlscoslty  model  hypothesizes  that  the 
Reynolds  stresses  are  linearly  related  to  the  mean  rate  of  strain.  The  rela¬ 
tionship  is  similar  to  that  derived  for  the  viscous  stresses: 


1  j  o  ij  2  ij, 

V  V  =  2v^e  -  "Y  g  k 


(5-1) 


where  denotes  the  Isotropic  eddy  viscosity,  which  depends  upon  the  local 
flow  conditions.  The  turbulent  kinetic  energy  term  has  been  added  to  make  the 
equation  valid  when  contracted  (multiplied  by  gij)^ 

-RljV^v^  =  2Vj.e^  -  2k  “  -  2k 


since  e^  =  0  for  an  Incompressible  fluid.  By  replacing  the  contravarlant 
strain  rate  with 

1  j  H 1  kj 

e  -  g  g 

*  I  +  V^,k> 

further  simplifications  yield 


(5-2) 


If  an  algebraically-specified  eddy-vlscoslty  model  Is  used,  substitution  of 
the  above  equation  Into  the  momentum  equation  (4-6)  completes  the  closure 
problem. 

In  the  k-e  model,  the  eddy  viscosity  is  related  to  the  turbulent  kinetic 
energy  and  Its  rate  of  dissipation  by 


since  k  and  e  are  scalars  which  are  Invariant  with  respect  to  coordinate 
transformation,  each  term  In  the  curvilinear-tensor  form  of  their  transport 
equations  has  a  direct  counterpart  In  the  Cartesian-tensor  form  of  these 
equations.  This  allows  one  to  simply  transform  each  term  of  the  modeled 
Cartesian  tensor  equations  to  the  curvilinear  tensor  form.  It  Is  assumed  that 
the  model  constants  are  not  affected  by  the  coordinate  transformation. 

The  k-G  model  equations  in  Cartesian  tensor  form  (see,  for  example, 
Bradshaw,  Cebecl,  and  Whltelaw  (1981))  are: 
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where  the  generation  term  Is 
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These  may  be  directly  converted  to  the  curvilinear  tensor  form: 
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(5-9) 

Note  that  the  convection  and  diffusion  terms  in  these  equations  have  the  same 
general  form  as  the  energy  equation  (3-45),  and  the  generation  term  G  is 
similar  to  the  dissipation  function  (3-44). 

Introducing  the  eddy-vlscoslty  relation  (5-2)  into  the  Reynolds  equations 
(4-9),  and  Including  the  continuity  equation,  gives  the  following  system  of 
six  coupled,  nonlinear  partial  differential  equations: 
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where  the  generally  accepted  model  constants  are:  “  0.09,  “  1.44, 

C  *  1.92.  a  ■  1.0,  and  a  -  1.3  .  These  equations  are  valid  for  Incom- 
e2  k  e 

presslble  turbulent  flow  In  any  coordinate  system.  No  approximations  have 
been  made  except  for  those  Inherent  In  the  turbulence  model.  Together  with 
the  proper  boundary  conditions,  equations  (5-10)  through  (5-13)  constitute  a 
well-closed  problem  which.  In  theory,  may  be  solved  for  the  six  unknowns;  V^, 
p ,  k ,  and  e . 


Equations  (5-10)  through  (5-13)  can  also  be  expanded  In  terms  of  their 
physical  components  to  yield: 
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VI.  REDDCTION  OF  BQDATIONS  TO  ORTHOGOHAL  COORDIMATES 

One  way  to  Illustrate  the  generality  of  the  equations  in  curvilinear 
tensor  form  is  to  simplify  them  to  the  more  familiar  equations  in  general 
orthogonal  coordinates  given  by  Nash  and  Patel  (1972).  The  physical 
components  of  velocity  must  be  used  when  attempting  to  compare  equations 
derived  by  tensor  methods  with  those  developed  using  orthogonal  curvilinear 
coordinates.  This  is  clearly  evident  because  the  equations  given  by  Nash  and 
Patel  are  written  using  the  physical  velocity  components.  This  exercise  also 
provides  an  example  of  the  use  of  tensor  analysis  and  the  introduction  of 
physically  measurable  vector  components  into  the  tensor  equations. 


VI. 1.  Slsplified  Tensor  Operations 

It  seems  useful  at  this  point  to  reiterate  several  of  the  tensor 
relations  which  take  simpler  forms  when  the  coordinate  system  is  orthogonal. 
Of  prime  Importance  is  the  metric  tensor  and  Its  Inverse: 


11  1 

®  .  2  (6-1) 

gjj  =  g^-^  =  0  for  i  j 


The  determinant  of  the  covariant  metric  tensor  or  Jacobian  Is 

g  =  =  (h^hjh^)^ 


Recall  that  (2-56)  gave  the  Chrlstoffel  symbols  for  an  orthogonal  coordinate 
system  as 


where  l,j,k  are  in  cyclic  order.  The  physical  components  of  a  tensor,  for 
Instance  the  contravarlant  velocity  vector,  can  be  found  from 

V(l)  =  (6-2) 


Nash  and  Patel  found  It  useful  to  introduce  the  curvature  parameter 
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which,  for  l^J  represents  the  curvature  of  the  surface  j=constant  In  the 
direction  1.  can  be  related  to  the  physical  components  of  the  Chrlstoffel 
symbols  In  the  following  manner; 
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where  l,J,k  are  In  cyclic  order.  The  other  geometric  coefficients  can  also  be 
simplified  by  noting  that 


g(ll)  »  1  (6-5a) 

g(lj)  =  0  (6-5b) 


Therefore,  D(im),  E(lm),  P(l),  Q(ln)  etc.  can  be  expressed  in  terras  of  the 
Chrlstoffel  symbols  or  curvature  parameters  as 
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+  (K^j^  +  Kjj^)  V(k)  =  0 


(6-15) 


where  l,j,k  are  In  cyclic  order.  The  reader  may  verify  that  this  result  Is 
Identical  to  the  continuity  equation  (2.11)  given  by  Nash  and  Patel  with 
(l,j,k)  =  (1,2,3)  and  V(l)  =  U,V,W. 


VI. 3.  Hoaentua  (Reynolds)  Equations 

If  we  Introduce  the  geometric  coefficients  of  (6-1)  through  (6-13)  Into 


the  Reynolds  equations  (4-8)  and  assume  F(l)  =  0,  the  mean  momentum  equations 
for  turbulent  flow  can  be  written  as 
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with  1,1, k  In  cyclic  order.  Letting  (5^,  5^,  5^)  =  (x,y,z)  and  (V(l),  V(2) , 


V(3))  =  (U,V,W),  and  then  rotating  the  Indices  1,1, k,  equations  (6-16)  become 
the  Reynolds  equations  (2.20)  through  (2.22)  In  Nash  and  Patel. 


VI. 4.  Reynolds  Stress  Transport  Equations 

The  transport  equation  (4-10)  for  the  Reynolds  stresses  v(l)v(1 )  may  be 
reduced  to  a  form  valid  for  an  orthogonal  coordinate  system  In  the  same  manner 
as  was  done  for  the  Reynolds  equations.  The  result  Is: 
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These  equations  are  Identical  to  equations  (2«28)  through  (2.30)  for  uv  (I*!, 
j=2,  k=3),  uw  (l=«3,j-l,k-2)  and  w  (1»2, j-3,k»l)  given  by  Nash  and  Patel  If 
the  equation  of  continuity  Is  employed  to  rearrange  the  triple  velocity  cor¬ 


relations. 


The  transport  equation  for  y  vCTTvCD  In  a  general  orthogonal  system  can 
be  derived  In  a  similar  manner  by  letting  l*j  In  equations  (4-10): 


l|  v(l)v(l)]  +  -5^^  v(l)v(l)]  +  Ij  v{l)v(l)] 

"l  ^  "j 

+  Ij  v(l)v(l)]  +  [KjjV(l)-K.^V(J)]  v(l)v(J)+[Kj^V(l)-K^jV(k)]  v(l)v(k) 

Ic  d  C 

+  v(l)v(l)  p  +  K.  .V(j)  +  K..  V(k)]  +  v(l)v(j)  p  -  K,,V(j)] 

"l  3C  J  j  35J 

+v(l)v(k)  -Kk<V(k)]+ 4  v(l)v(l)v(l)l+  (1  v(l)v(l)v(j)] 

k  35  1  35^  ^  "j  3CJ  ^ 


1  3  ,1 


[-0  v(l)v(l)v(k)]  +  TT  (K..+  K,  .)  v(l)v(l)v(l)  +  (■»  K,. 


+  K,  .)  v(l)v(l)v(1)  +  (-T  K.,  +  K.,  )  v(l)v(l)v(k) 


v(l)v(j)v(j) 


K.  .v(l)v(k)v(k)  + 


.  v(^)  3v(J)  v(l)  3v(.1)  »  v(i)  3v(k 

\  "j  3cJ  “‘"1  35* 


v(l)  3v(k 

\  35" 


+  v(l)v(l)  +  v(l)v(J )  +  v(l)v^k)}  =  0 


(6-18) 


When  the  continuity  equation  Is  employed,  equation  (6-18)  reduces  to  equations 
(2.24)  through  (2-26)  In  Nash  and  Patel  for  u^,  v^,  and  w^. 


VII.  EVALDATIOH  OF  TYPICAL  CHRISTOFFEL  SYMBOLS  FOR 
ROHERICALLY-GEHERATED  COORDINATES 


For  most  practical  applications  requiring  the  use  of  curvilinear  tensor 
equations,  it  is  most  unlikely  that  an  analytical  coordinate  system  will  be 
available.  The  required  geometric  coefficients  must  usually  be  evaluated  from 
numerically-generated  coordinates.  Naturally  the  question  Is  asked  as  to  what 
sort  of  geometric  coefficients  are  produced  from  a  typical  grid  generation 
method.  To  Illustrate  the  answer  a  few  typical  Chrlstoffel  symbols  for  two, 
shlp-llke  three  dimensional  bodies  will  be  shown.  The  magnitude  and  variation 
of  the  Chrlstoffel  symbols  are  of  interest  because  of  their  ubiquity  in  the 
equations  of  fluid  motion. 

The  two  coordinate  systems  considered  are  those  for  the  3:1  elliptical 
cross-section  body  of  Groves  et  al.  (1982)  and  the  SSPA  cargo  ship  hull  of 
Larsson  (1974).  Flow  computations  for  these  bodies  using  partlally-tans- 
formed,  partially-parabollc  equations  In  numerically-generated,  body-fitted 
coordinates,  have  been  carried  out  by  Chen  and  Patel  (1985a,b).  The 
Chrlstoffel  symbols  shown  herein  were  calculated  from  the  numerical  coordin¬ 
ates  used  by  Chen  and  Patel  In  their  computations  and  may  be  considered  as 
representative  of  the  geometric  coefficients  which  may  be  encountered  in 
practice. 

The  numerical  coordinate  grids  for  these  bodies  were  generated  by  using 
the  Poisson  equation  technique  described  by  Chen  and  Patel  (1985a).  The 
coordinate  system  and  general  orientation  of  the  base  vectors  are  shown  In 
figure  3.  Different  views  of  the  coordinate  grids  for  the  3:1  body  and  the 
SSPA  hull  are  shown  In  figures  4  and  5,  respectively. 

The  relationship  between  the  Chrlstoffel  symbols  and  the  base  vectors  can 
be  readily  Identified  by  comparing  (2-52)  and  (2-58),  from  which  one  concludes 
that 


=  ^jk 


(7-1) 


Hence,  the  Chrlstoffel  symbol  represents  the  rate  of  change  of  the  jth  base 
vector  with  respect  to  the  kth  coordinate  resolved  along  the  1th  coordinate 
direction.  The  Chrlstoffel  symbols  can  be  calculated  from  (2-54): 
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Evidently,  the  Chrlstoffel  symbols  depend  upon  the  metric  coefficients  and 


and  vere  calculated  using  standard 


their  derivatives.  The  gj[j ,  » 

finite-difference  techniques.  Derivatives  were  evaluated  using  central  dif¬ 
ferences  wherever  possible  and  with  forward  or  backward  differences  on  the 
domain  boundaries. 


2 

Figures  6a  and  6b  show  the  variation  of  on  the  n  =  2,3  coordinate 

surfaces  of  the  3:1  body.  Forward  differencing  seems  to  be  the  cause  of  the 

2 

jagged  variation  on  n  =  2  .  The  variation  of  is  much  smoother  on  the  next 

coordinate  surface,  n  *  3,  where  central  differencing  can  be  used  to  calculate 

the  necessary  derivatives.  Recall  that,  for  an  orthogonal  coordinate  system, 
2 

r._  is  zero.  The  Influence  of  the  rapid  variation  of  the  body  geometry  is 

3  2 

shown  in  figures  6c  through  6h.  Near  the  body  and  show  more  variation 
compared  to  the  regions  well  removed  from  the  body  surface.  This  is  particu¬ 
larly  evident  at  C  =  40,  where  the  body  terminates  and  the  coordinate  lines 
make  a  transition  from  elliptical  to  nearly  circular  shapes  (see  figures  4d 
and  4e).  In  the  wake  region  the  Chrlstoffel  symbols  again  have  a  smooth 
variation.  Thus,  whenever  the  geometry  of  the  bounding  surface  changes  rapid¬ 
ly  one  must  expect  equally  rapid  changes  in  the  geometric  coeffclents. 


Further  effects  of  the  body  geometry  on  the  Chrlstoffel  symbols  can  be 

3 

seen  in  figure  7  which  shows  the  variation  of  r  _  for  the  SSPA  hull.  Note 

^  3 

that  the  hull  terminates  at  C  =  20.  The  variation  of  in  this  case  is 

qualitatively  similar  to  that  for  the  3:1  body  except  near  the  hull  surface. 

Because  flow  calculations  using  the  complete  curvilinear  tensor  formula¬ 
tion  of  the  governing  equations  have  only  just  begun,  it  is  not  possible  at 
this  time  to  fully  assess  the  Influence  of  the  magnitude  and  variation  of  the 
geometrical  coefficients  upon  the  calculations.  If  the  so-called  wall  func¬ 
tions  are  used  to  satisfy  the  boundary  conditions  on  the  body  for  turbulent 
flows,  the  coordinate  surface  n  -  2  is  not  used  in  the  calculations.  Clearly 
this  alleviates  some  difficulties  because  the  coefficients  on  n  =  3  can  be 
evaluated  using  central  differencing  and  figure  6b  shows  that  F^^  has  a  smooth 
distribution.  Trouble  may  arise  if  the  equations  must  be  Integrated  to  the 
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wall,  as  will  be  the  case  for  laminar  flows,  or  If  low-Reynolds-number  turbu¬ 
lence  models  are  employed.  At  the  outer  boundaries  the  problems  may  not  be  as 
severe.  One  can  simply  specify  the  boundary  conditions  'one  step'  Inside  the 
domain  boundary,  unless  of  course  computations  are  being  done  for  a  duct 
flow.  Observe  that  the  Chrlstoffel  symbols  are  not  the  only  geometric  coeffi¬ 
cients  which  show  up  In  the  equations.  Combinations  of  the  metric  tensor, 
Chrlstoffel  symbols,  and  derivatives  of  the  Chrlstoffel  symbols  appear 
throughout  the  equations.  Obviously  additional  study  Is  required  to 
completely  explore  the  Influence  of  the  geometric  coefficients  upon  the 
quality  of  flow  computations. 


APPENDIX 

Vector  Operatloos  In  Teosor  Form 


The  usual  vector  operations  can  be  readily  expressed  In  terms  of  curvi¬ 
linear  tensors  using  the  methods  outlined  In  Section  II.  The  results  are 
summarized  below. 


The  scalar  product  of  two  vectors  is 


or  In  terms  of  physical  components. 


g^^A(l)B(]) 


COS  A(l)  B(j) 


(A-1) 


(A-2) 


The  vector  product  which  yields  a  contravarlant  vector  Is 

;ljk 
“ 

g 

or,  defining  a  covariant  vector, 


e  1 

A.B,* 


A  X  g  =  -J72 


(A-3) 


A  X  5  -  Cj 


(A-4) 


The  reader  may  verify  that  this  reduces  to 

C(l)  -  Cijj^  A(J)B(k)  (A-5) 

for  an  orthogonal  coordinate  system. 

The  gradient  of  a  scalar  function  forms  the  components  of  the  covariant 
vector  given  by 


By  introducing  the  physical  components  this  becomes 


Vt(l)  ”  (no  sum  on  1) 


for  general  coordinates,  and  reduces  to 


V^(l) 


1 

r - ^  (no  sum  on  1) 

1  3C 


for  orthogonal  coordinates. 


The  divergence  of  a  vector  la  obtained  by  contracting  the 
derivative  of  a  contravarlant  vector, 

dlv  (A)  =•  V  •  ^  ^  ^Ik 

,1 

or,  using  (2-55), 


V«A 


1  3 

IT72  .-1 


A^) 


The  curl  of  a  vector  can  he  expressed  by  either 


V  X  A 


.ijk 

772 


or 


V  X  A 


.IJk 

T72  ®jp 


The  Laplaclan  of  a  scalar  function  (j> ,  defined  as 

2 

V  ^  =  V*  (ViJ)  ) 


(A-7) 


(A-8) 


covariant 


(A-9) 


(A-10) 


(A-lla) 


(A- 11b) 


(A-12) 


Is  a  scalar  given  by 


„2,  mn, 

7  <|)  -  g  (fi 


^«n.mKnn>  afW 

3C  o€ 


1/2  mn 


g(mn)  -  g 


f(in)  -  -  g(rs)r(inrs) 


a  _  -1/2  a 

ac (m)  "  ®mm  m 
3s 


(A-13) 


(A-14) 


(A- 14a) 


(A-14b) 


(A-14c) 


Substituting  the  physical  components  for  the  tensor  components  yields  the 
following  Laplaclan  of  a  vector; 

2  mn  1  -1/2  ^^2  .  .  ,  .  aA(l)  ,  .  aA(s) 

^  ^  =  g^^  {7  A(l)  +  2E(lm)  +  2C(lms) 


+  P(1)A(1)  +  Q(ls)A(s)} 


where 


E(lm)  =  g(mn)D(ln) 


C(lms)  =  g{mn)r(lns) 


(A-15) 


(A- 15a) 


(A- 15b) 


P(i)  =  (A-15c) 

Q(ls)  =  ^ 3 ^ (m)'^  *  D(8m)  +  D(mm)  +  r(nnm)]  C(lms) 


+  C(lmn)  r(nras) 


D(lm) 


1/2  ^^11  ^1 _  ^^11 

®11  ac(m)  “  2g^j  ac(ra) 


(A-15d) 


(A-15e) 


■  .  *  .  ■  .  •  **■»“*.' 


yj'js 


«■  V  V ' 


i 


kvS> 

r  -  *  - 


Green's  theorem  written  In  curvilinear  tensor  form  is 


/y  dV  -  /g  A^n^dS  (A-16) 

where  n^^  is  the  outward  unit  normal  to  the  surface  S  bounding  the  volume  V. 
Other  forms  of  Green's  theorem  are 

°  U  ^1  “  -^S 

where  ^  is  the  intrinsic  derivative  normal  to  the  surface  S. 

Stokes*  theorem  is  given  by 
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Figure  5.  Continued 

(b)  Keel  View. 
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Figure  5.  Continued 

(e)  Transverse  View  at  x/L  =  2.302. 
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